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Abstract 

The construction of specific supersymmetric grand unified models based on the 



00 

o 

j-». ' Pati-Salam gauge group and leading to a set of Yukawa quasi-unification conditions is 

^^ , briefly reviewed. For each sign of the parameter ^, an appropriately chosen condition 

c7^ ' from this set can allow an acceptable 6-quark mass within the constrained minimal 

^^ , supersymmetric standard model. The restrictions on the parameter space which arise 

O ' from the cold dark matter constraint, the inclusive decay b — > s'y, the muon anomalous 

magnetic moment and the collider bounds are also investigated. For /x > 0, a wide 

and natural range of parameters is allowed. On the contrary, the ^ < case not 

(~| , only is disfavored by the present experimental data on the muon anomalous magnetic 

moment, but also is excluded by the combination of the cold dark matter and 6 ^ 57 

requirements. In the /i > case, the predicted neutralinos are possibly detectable in 

^\. ' the future direct cold dark matter searches. Moreover, the fj, term is generated via a 

Peccei-Quinn symmetry and proton is practically stable. The same model gives rise to 
a new shifted inflationary scenario, which is based only on renormalizable terms, does 
not suffer from the problem of monopole overproduction at the end of inflation and is 
compatible with the cosmic microwave background radiation constraint. Although the 
relevant part of inflation takes place at values of the inflaton fleld which are not much 
smaller than the reduced Planck scale and, thus, supergravity corrections could easily 
invalidate inflation, this does not happen provided that an extra gauge singlet with a 
superheavy vacuum expectation value, which originates from D-terms, is introduced 
and a speciflc form of the Kahler potential is used. The constraint from the cosmic 
background explorer can again be met by readjusting the values of the parameters 
which were obtained with global supersymmetry. 
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1 Introduction 

The constrained minimal supersymmetric standard model (CMSSM), which is a 
highly predictive version of the minimal supersymmetric standard model (MSSM) 
based on universal boundary conditions [1], can be further restricted by being em- 
bedded in a supersymmetric (SUSY) grand unified theory (GUT) with a gauge group 
containing S'C/(4)c and SU{2)ji. This can lead [2] to 'asymptotic' Yukawa unifica- 
tion (YU) [3], i.e. the exact unification of the third generation Yukawa coupling 
constants at the SUSY GUT scale Mgut- Indeed, assuming that the electroweak 
Higgs superfields Hi, H2 and the third family right handed quark superfields f^, 
6^ form SU{2)ii doublets, we obtain [2] the asymptotic Yukawa coupling relation 
ht = hb and, hence, large tan/3 ~ mt/nib. Moreover, if the third generation quark 
and lepton SU{2)l doublets [singlets] q^ and ^3 [b'^ and t'^] form a S'[/(4)c 4-plet 
[4-plet] and the Higgs doublet Hi which couples to them is a 5[/(4)c singlet, we 
get hb = hr and the asymptotic relation nib = itt-t follows. The simplest GUT 
gauge group which contains both 5C/(4)c and SU{2)ji is the Pati-Salam (PS) group 
Gps = SU{4)c X SU{2)l X SU{2)r (for YU within SO{W), see Ref. [4]). 

However, given the experimental values of the top-quark and tau-lepton masses 
(which, combined with YU, naturally restrict tan/? ~ 50), the CMSSM supple- 
mented by the assumption of YU yields unacceptable values of the 6-quark mass 
for both signs of the parameter fi. This is due to the generation of sizeable SUSY 
corrections [5] to nib (about 20%), which arise from sbottom-gluino (mainly) and 
stop-chargino loops [5, 6, 7] and have the sign of fi (with the standard sign conven- 
tion of Ref. [8]). The predicted tree-level nib{Mz), which turns out to be close to 
the upper edge of its 95% confidence level (c.l.) experimental range 

2.684 GeV < mb{Mz) < 3.092 GeV (with as{Mz) = 0.1185), (1) 

receives, for // > [/U < 0], large positive [negative] corrections which drive it well 
above [a little below] the allowed range. This range is derived [9] from the 95% c.l. 
range for mb{nib) in the MS scheme (3.95 — 4.55 GeV) [10] evolved up to Mz in 
accord with the analysis of Ref. [11]. Consequently, for both signs of fi, YU leads 
to an unacceptable mb{Mz)-, with the ^ < case being much less disfavored. 

The usual strategy to resolve this discrepancy is the introduction of several kinds 
of nonuniversalities in the scalar [12, 13, 14] and/or gaugino [15, 16] sector of MSSM 
with an approximate preservation of YU. On the contrary, in Ref. [9], this problem 
is addressed in the context of the PS GUT, without need of invoking departure from 
the CMSSM universality. The Higgs sector of the model is extended by including 
an extra SU{4)c non-singlet Higgs superfield with Yukawa couplings to the quarks 
and leptons. The Higgs SU{2)l doublets contained in this superfield can naturally 
develop [17] subdominant vacuum expectation values (VEVs) and mix with the main 
electroweak doublets which are assumed to be S'C/(4)c singlets and form a SU{2)r 
doublet. This mixing can, in general, violate SU{2)r. Consequently, the resulting 
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electroweak doublets Hi, H2 do not form a SU{2)fi doublet and break SU{4:)c too. 
As a consequence, a moderate violation of the YU is obtained, which can allow 
an acceptable 6-quark mass even with universal boundary conditions. Obviously, a 
small deviation from YU is enough for an acceptable prediction of the 6-quark mass 
when /U < 0, while, for ^ > 0, a more pronounced deviation is needed. 

In this review, we outline the construction of two particular SUSY GUT models 
which can cause a deviation from YU adequate for /i > [/i < 0] in the first 
[second] model. We then discuss the resulting CMSSM in each case and the various 
astrophysical and experimental requirements which restrict its parameter space. 
They originate from the data on the cold dark matter (CDM) abundance in the 
universe derived by the Wilkinson microwave anisotropy probe (WMAP) [18, 19], 
the inclusive branching ratio of 6 ^ 57 [20], BR(6 -^ 57), the muon anomalous 
magnetic moment a^ [21] and the mass of the lightest Higgs boson nih [22]. We 
show that, for /i > 0, our model possesses a wide range of parameters which is 
consistent with all these constraints. On the contrary, for fi < 0, the upper bound 
on the mass of the lightest sparticle (LSP) from the CDM abundance in the universe 
is incompatible with the data on BR(6 -^ 57). Thus, the latter scheme is not viable. 

The /i > model possesses a number of other interesting features too. The LSP is 
possibly detectable in the near future direct CDM searches. The fi problem of MSSM 
is solved [23] via a Peccei-Quinn (PQ) symmetry [24] which also solves the strong CP 
problem. Although the baryon (i?) and lepton (L) numbers are explicitly violated, 
the proton lifetime is considerably higher than the present experimental limits. Fur- 
thermore, the same model leads to a new shifted hybrid inflationary scenario [25], 
which avoids monopole overproduction at the end of inflation and reproduces the 
results on the quadrupole anisotropy of the cosmic microwave background radiation 
(CMBR) from the cosmic background explorer (COBE) measurements [26]. The 
relevant part of inflation takes place at values of the inflaton field which are not 
much smaller than the Planck scale and, thus, supergravity (SUGRA) corrections 
could easily invalidate it. It is, however, shown that infiation can be kept intact 
provided that an extra gauge singlet is introduced and a specific form of the Kahler 
potential is used. Although the SUGRA corrections are sizable, the constraints from 
COBE can again be met by readjusting the values of the parameters which were 
obtained with global SUSY. 

The construction of the models is briefly reviewed in Sec. 2 and the resulting 
CMSSM is presented in Sec. 3. The parameter space of the CMSSM is restricted 
in Sec. 5 taking into account a number of cosmological and phenomenological re- 
quirements which are exhibited in Sec. 4. The deviation from YU is estimated in 
Sec. 6. Issues related to the direct LSP detection are examinated in Sec. 7. The 
resolution of the fi problem of MSSM and the stability of the proton are discussed 
in Sees. 8 and 9 respectively. The new shifted inflationary scenario is outlined in 
Sec. 10. Finally, we summarize our conclusions in Sec. 11. Details concerning the 
evaluation of the neutralino— nucleus cross section are given in Appendix A. 
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2 The Pati-Salam SUSY GUT Models 



2.1 The General Set-up 

We focus on a SUSY GUT model based on the PS gauge group Gps = 5C/(4)c x 
SU{2)l X SU{2)r described in detail in Ref. [27] (see also Ref. [28]). The repre- 
sentations and the transformation properties under Gps as well as the extra global 
charges of the various superfields contained in the model are presented in Table 
1, which also contains the extra Higgs superfields required for accommodating an 
adequate violation of YU (see below). 

The left handed quark and lepton superfields of the rth generation (r = 1, 2, 3) 
are accommodated in a pair of superfields 
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where tilde denotes transpose with respect to (w.r.t.) SU{A)c- The gauge symmetry 
Gps can be spontaneously broken down to the standard model (SM) gauge group 
(Csm) through the VEVs which the superfields 
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(3) 



acquire in their right handed neutrino directions v'fj and U'fj. The model also contains 
a gauge singlet S which triggers the breaking of Gps, a 5f/(4)c 6-plet G which gives 
[29] masses to the right handed down quark type components of H'^, H^ and a pair of 
gauge singlets N , N for solving [23] the /x problem of the MSSM via a PQ symmetry 
(see Sec. 8). In addition to Gps, the model possesses two global U{1) symmetries, 
namely a PQ and a R symmetry, as well as a discrete Z'^'^ symmetry ('matter parity') 
under which F, F'^ change sign. Note that global continuous symmetries such as 
our PQ and R symmetry can effectively arise [30] from the rich discrete symmetry 
groups encountered in many compactified string theories (see e.g. Ref. [31]). 

In the simplest realization of this model [29] , the electroweak doublets Hi , if 2 
are exclusively contained in the bidoublet superfield /i, which can be written as 



h 



(^^"OM^i :n 



(4) 



Under these circumstances, Hi = hi with i = 1,2 and so the model predicts YU 
at Mqut (which is determined by the requirement of the unification of the gauge 
coupling constants), i.e. 

ht:h^:hr = 1:1: I, (5) 
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Table I : The representations and transformations under Gps as well as the extra global charges 
of the superfields of our model (Uc G S'[/(4)c, Ul G SU{2)l, Ur G SU(2)r and T ,t and * stand 
for the transpose, the hermitian conjugate and the complex conjugate of a matrix respectively). 



since, after the breaking of Gps to GsM, the third family fermion masses originate 
from a unique term of the underlying GUT as follows: 



y33F3{h)Fi 



2/33 



-V2 tf + VihJf + ViT t") + 



where Vi = {Hi). (6) 



A moderate violation of YU can be accommodated by adding two new Higgs 
superfields h' and h! with 



h 



' = [ K K ) and h' = f h'2 h[ ] 



(7) 



In accordance with the global symmetries imposed (see Table 1), h' can couple 
to FF'^ since FF'^ = (15,2,2), whereas h' can give mass to the color non-singlet 
components of h' through a superpotential term mh'h' {m ~ Mqut — 2x 10"*^^ GeV), 
which corresponds to the following Gps invariant quantity: 



mir 



(h'eh'^e) 
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and Tr denotes trace taken w.r.t. the S'C/(4)c and SU{2)l indices. Other important 
superpotential terms, which, being non-renormahzable, are suppressed by the string 
scale Ms ~ 5 X 10^"^ GeV, are 



smce 



iJ'=i?^/i'/i= 1, 1, 1 X 1 + 3 X 3 , (9) 



H'H' = (4, 1, 2)(4, 1, 2) = (15, 1, 1 + 3) + • • • , 
/I'/i = (15, 2, 2)(1, 2, 2) = (15, 1, 1 + 3) + • • • . 

We see that, in Eq. (9), there are two independent couphngs: a couphng between 
the SU{2)ji singlets in H'^H^ and h'h, {H^H^)\h'h, and a coupling between their 
triplets, (H'^H^)sh'h. As it turns out, the singlet coupling provides us with an 
adequate deviation from YU for /i < 0. On the other hand, ^ > requires a 
stronger deviation. 

2.2 The /x > Case 

The necessary deviation from YU for ^ > can be obtained by a further enlargement 
of the Higgs sector so that contributions to the coupling between h and h' from 
renormalizable terms are also allowed. To this end, we introduce two additional 
Higgs superfields (p, (p. The superfield (j) is aimed to give superheavy masses to the 
color non-singlets in through a term (j)(j), whose coefficient is of order Mqut- The 
superfield (p, on the other hand, yields the unsuppressed coupling X^cph'h (with A3 
being a dimensionless constant), which overshadows the coupling {H'^H'^)^h'h and 
corresponds to the Gps invariant term 

AsTr (h'ecph^e) . (10) 



During the spontaneous breaking of Gps to GsM, <t> acquires VEV in the SM 
singlet direction. Therefore, 

(0) = z;^ (^l^ 1, -^) , (11) 

where v^ ~ Mqut and the structure of {(j)) w.r.t. Gps is shown with 

Ti5 = -i^diag(l,l,l,-3) and (73 = diag(l, -l) . (12) 

Expanding the superfields in Eq. (7) as linear combination of the fifteen genera- 
tors T" of S'[/(4)c with the normalization Tr(T"T^) = (5"^ and denoting the SM 
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singlet components with the superfield symbol, we can easily establish the following 
identities: 



Tr (h'eh'^ej = h'^eh'^ + h'^eh'^ + ■■■, (13) 

Tr (h'e {T^\ 1, ag) h" e) = (/ife/is - hjeh'^ , (14) 



where the notation of Eq. (11) has been applied. 

Inserting Eq. (11) in Eq. (10), employing Eqs. (13) and (14), and collecting 
Eqs. (8) and (10) together, we get the mass terms 

mh'^e {h'2 + 02/12) + rn (h'^ + ai/i]") eh'2, (15) 

where the mixing effects are included in the following coefficients: 

"1 = -012 = -\sv,j,/\/2m. (16) 

2.3 The ^Ji <Q Case 

In the n < case, an adequate violation of YU can be achieved (without the inclu- 
sion of (p and (p) predominantly via the non-renormalizable SU{2)pi singlet coupling 
X-i_{H'^H'^)ih'h/Ms (with Ai being a dimensionless constant), which corresponds to 
the Gps invariant term 

^Tr (h'e{H^^H^^We) . (17) 

During the spontaneous breaking of Gps to Gsm, H'^ and H'^ acquire VEVs in the 
SM singlet direction. Therefore, 



{H^^){H^^))^ = vi,.{THc,l,^), (18) 

where vh<^ ~ AfcuT and the notation of Eq. (11) has been applied with 

Tffc =diag('o,0,0, 1) and do = diagf 1, iV (19) 



Expanding the superfields in Eq. (7) as we did in deriving Eq. (14), we obtain the 
following identity: 

Tr (h'e {Thc, 1, ctq) h^e) = -Vs (h'^eh2 + hjeh'^) /2. (20) 



Inserting Eq. (18) in Eq. (17), employing Eqs. (13) and (20), and collecting Eqs. (8) 
and (17) together, we end up again with Eq. (15) but with 

which are suppressed by Mgut/-^s- 
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2.4 Yukawa Quasi-Unification Conditions 

It is obvious from Eq. (15) that we obtain two combinations of superheavy massive 
fields 



h[, H[ and h'2, H2, where H'^ = — f , , = , i 



^/T+\a■ 



.12 • 



1,2. 



(22) 



The electroweak doublets Hi, which remain massless at the GUT scale, are orthog- 
onal to the H' directions: 



Hi 



-a*h'i + hi 



^/l + 



ai 



Solving Eqs. (22) and (23) w.r.t. hi and h[, we obtain 

-a^H, + m 



hi 



^ + «*^ and/.' 



Vi + RF 



Vi + RF 



(23) 



(24) 



Consequently, the third family fermion masses are now generated by the following 
terms (compare with Eq. (6)): 



y33F3{h)Fi + 2y'^^Fs{h')R 



I l-pa2/V3 ,„ l — pai/y/3 , ,^ 1 + VSpai A 



\[^ 



"2 



vi + KF 



\/l + |ai| 



where p = ^33/^33 with < p < 1 and the color non-singlet components of h'^ 
are included in the ellipsis. The fields H[, being superheavy, contribute neither to 
the running of renormalization group equations (RGEs) nor to the masses of the 
fermions. Consequently, the asymptotic exact YU in Eq. (5) can be replaced by a 
set of asymptotic Yukawa quasi-unification conditions (YQUCs) : 



ht : hb : K 



1 - /9Q2/\/3 



Vl + I«2p 



1 — pa\l\f?) 



vt+kF 



1 + \f?>pa\ 



vt+kF 



(25) 



Substituting Eqs. (16) and (21) into Eq. (25), we obtain 

(1 + c) : (1 - c) : (1 + 3c) with < c < 1 for ^ > 



^*-^^-^""|(l-c):(l-c) :(l + 3c) with -1/3 < c < for ^ < ' ^^^^ 

For simplicity, we restricted ourselves to real values of c only, where c = pai/^/?,. 

The deviation from YU can be estimated by defining the following relative split- 
tings: 

hb[T] -^t _i -2c/ {I + c) [2c/(l + c)] for p > 



5h 



6[r] 



ht 



[Ac/ {I + c)] 



for p < 



(27) 
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Figure 1: A running of the Yukawa coupling constants from Q = Mgut toQ — Mz for nib(Mz) 
2.888 GeV, p > 0, ttilsp = 200 GeV and Af^ = 1 corresponding to tan/3 ~ 58. 
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Figure 2: A ruruiing of the Yukawa coupling constants from Q — Mgut to Q — Mz for mb{Mz) 
2.888 GcV, ^ < 0, rriLSP = 350 GeV and Af^ = corresponding to tan/3 ~ 47. 
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3 The resulting CMSSM 

Below MguT) the particle content of our models reduces to this of MSSM (modulo 
SM singlets). We assume universal soft SUSY breaking scalar masses rriQ, gaugino 
masses M1/2 and trilinear scalar couplings ^0 at Mqut- Therefore, the resulting 
MSSM is the so-called CMSSM [1] supplemented by a suitable YQUC from the set 
in Eq. (26). Let us emphasize that the specific relations between qi and 02 given 
in Eq. (16) [Eq. (21)] for ^u > [/U < 0] ensure a SUSY spectrum which leads to 
successful radiative electroweak symmetry breaking (REWSB) and a neutral LSP 
in a large fraction of the parametric space in the framework of the CMSSM (and 
not in a general version of MSSM). 

We integrate the two-loop RGEs for the gauge and Yukawa coupling constants 
and the one-loop ones for the soft SUSY breaking terms between Mqut and a 
common SUSY threshold MsusY — ("^ti^T-ta)^ (^1.2 are the stop mass eigenstates) 
determined in consistency with the SUSY spectrum. At MsusY, we impose REWSB, 
evaluate the SUSY spectrum and incorporate the SUSY corrections to the b and r 
mass [6, 7, 12]. The corrections to mr (almost 4%) lead [32] to a small decrease 
[increase] of tan f3 for /i > [;U < 0] . From MgusY to Mz , the running of gauge and 
Yukawa coupling constants is continued using the SM RGEs. 

For presentation purposes, M1/2 and rriQ can be replaced [32] by the LSP mass, 
rriLSP) and the relative mass splitting, A^j, between the LSP and the lightest stau 
f2. For simplicity, we restrict this presentation to the ^0 = case (for Aq ^ 0, see 
Refs. [9, 33]). So, our free input parameters are 

sign(/i), rriLSP, A^2, with Af^ = {mf^ - mLSp)/'mLSP- 

For any given mh{Mz) in the range in Eq. (1) with fixed masses for the top quark 
fn,t{mt) = 166 GeV and the tau lepton mr{Mz) = 1.746 GeV, we can determine the 
parameters c and tan (3 at MsusY so that the corresponding YQUC in Eq. (26) is 
satisfied. Consequently, for any given mh{Mz), a prediction for tan/3 can be made, 
in contrast to the original version of CMSSM [1]. 

In Fig. 1 [Fig. 2], we present a Yukawa coupling constant running from Mqut 
to Mz for the central value of mb{Mz) = 2.888 GeV. At Mqut, we apply Eq. (26) 
with c = 0.154 [c = —0.043] (corresponding to tan/3 ~ 58 [tan/3 ~ 47] at Msusy) 
for ^ > [/^ < 0], rriLSP = 200 GeV [tolsp = 350 GeV] and A^^ = 1 [A^^ = 0]. 
The kinks on the various curves correspond to the point where the MSSM RGEs 
are replaced by the SM ones. We observe that, for /j, > 0, hr and h{j split from ht by 
the same amount but in opposite directions with h^ becoming smaller than ht (see 
Fig. 1), while, for fi < 0, ht and hi, remain unified (see Fig. 2). 

In Fig. 3 [Fig. 4], we present the values of the mass parameters niA, ttiq, Mi 12 
and Msusy versus to-lsp for ^u > [/U < 0], A^j = 1 [Afj = 0] (see Sec. 5.1 [Sec. 5.2]) 
and with mh{Mz) = 2.888 GeV. We observe that tuq > M1/2 [mo < M1/2] and 
ruA < 2?n,LSP i^T^A ^ 2r7T,LSp] for ^u > [;U < 0] . The relation between m-LSP and uia 
or rrif^ is crucial for the reduction of the LSP relic abundance (see Sec. 4.1). 
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4 Cosmological and Phenomenological Constraints 

The parameter space of our models can be restricted by using a number of phe- 
nomenological and cosmological constraints. We will now briefly discuss these re- 
quirements (for similar recent analyses, see Refs. [14, 15, 16, 34]). 

4.1 Cold Dark Matter Considerations 

According to the WMAP results [18], the 95% c.l. range for the CDM abundance is 

f^CDM/i' = 0.1126+S. (28) 

In the context of the CMSSM, the LSP can be the lightest neutralino x- It naturally 
arises [35] as a CDM candidate. We require its relic abundance, Q.i.sph'^, not to 
exceed the upper bound derived from Eq. (28) (the lower bound is not considered 
since other production mechanisms of LSPs may be present too [36, 37] and/or other 
CDM candidates [38, 39] may also contribute to ^cDMh'^)'- 

f^LSp/i' < 0.13. (29) 

For each sign of /i, an upper bound on rriLSP (or m^) can be derived from Eq. (29). 
We calculate JIlsp^^ using the publicly available code micrOMEGAs [40] (not the 
latest version [41]). This includes accurately thermally averaged exact tree- level 
cross sections of all possible (co) annihilation processes, treats poles properly and 
uses one-loop QCD (not SUSY QCD [41, 43]) corrections to the Higgs decay widths 
and couplings to fermions [42]. Good agreement between this code and other inde- 
pendent calculations of Olsp^i^ which include the vl-pole effect and neutralino-stau 
coannihilation is recently reported in the first paper of Ref. [9] and in Ref. [43]. 

In most of the parameter space of the CMSSM, the LSP is an almost pure bino 
and r^LSP^^ increases with rriLSP- Therefore, Eq. (29) sets a very stringent upper 
limit on the LSP mass [44]. However, as pointed out in Refs. [45, 46], a substantial 
reduction of f^LSP^^ can be achieved in some regions of the parameter space thanks 
to two mechanisms: the A-pole effect and the bino-slepton coannihilations. The first 
is activated for tan /? > 40 [tan /? ~ 30 — 35] for ^u > [/U < 0] , where the presence of 
a resonance (2?n,LSP — ^tt-a) in the LSP annihilation via a s-channel exchange of an 
^-boson is possible. On the other hand, coannihilations can be activated for every 
tan P and both signs of fi, but this requires a proximity between the masses of the 
LSP and the next-to-LSP, which turns out to be the f2 for tan/3 > 10 [46] and not 
too large values of ^o [47] or ruQ [48]. For fixed rriLSP) ^Lsph'^ decreases with A^j, 
since the x ~ ^2 coannihilations become more efficient. So the CDM criterion can 
be used for restricting Afj . 

Our models give us the opportunity to discuss how both these reduction mech- 
anisms operate. As we noticed from Figs. 3 and 4, for ^ > 0, there is a significant 
region dominated by the ^-pole effect, while, for /U < 0, the x ~ ^2 coannihilation is 
the only available reduction mechanism. 
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4.2 Branching Ratio of 6 — > 57 

Taking into account the recent experimental results [20] on BR(6 -^ 57) and com- 
bining [9] appropriately the various experimental and theoretical errors involved, we 
obtain the following 95% c.l. range: 

a) 1.9 X 10"^ < BR(6 -^ 57) and b) BR(6 -^ 57) < 4.6 x 10"^. (30) 

We compute BR(6 -^ 57) by using an updated version of the relevant calculation 
contained in the micrOMEGAs package [40]. In this code, the SM contribution is calcu- 
lated using the formalism of Ref. [49]. The charged Higgs boson, H^, contribution 
is evaluated by including the next-to-leading order (NLO) QCD corrections from 
Ref. [50] and the tan/3 enhanced contributions from Ref. [51]. The dominant SUSY 
contribution, BR(6 -^ ■s7)|susY, includes resummed NLO SUSY QCD corrections 
from Ref. [51], which hold for large tan/3. The H contribution interferes construc- 
tively with the SM contribution, whereas BR(6 -^ >s7)|susY interferes destructively 
[constructively] with the other two contributions for fi > [/j, < 0]. Although the 
improvements of Ref. [52] are not included in this routine, their impact is not im- 
portant for yU > 0, whereas, for /i < 0, they are not expected to change essentially 
our conclusions. The SM plus H^ contribution and the BR(6 -^ s7)|susY decrease 
as ?n,LSP increases and so a lower bound on rriLSP can be derived from Eq. (30a) 
[Eq. (30b)] for ^ > [/U < 0] with the bound for ^u < being much more restrictive. 

4.3 Muon Anomalous Magnetic Moment 

The deviation, 5a^, of the measured value of the muon anomalous magnetic moment, 
a^, from its predicted value in the SM, a^, can be attributed to SUSY contributions 
arising from chargino-sneutrino and neutralino-smuon loops. The quantity 5a ^ is 
calculated by using the micrOMEGAs routine based on the formulas of Ref. [53]. 
The absolute value of the result decreases as rriLSP increases and its sign is positive 
[negative] for // > [/x < 0]. 

On the other hand, the calculation of afj^ is not yet stabilized mainly because 
of the instability of the hadronic vacuum polarization contribution. According to 
the most up-to-date evaluation of this contribution in Ref. [54], there is still a 
discrepancy between the findings based on the e^e^ annihilation data and the ones 
based on the r-decay data. Taking into account these results and the recently 
announced experimental measurements [21] on a^, we impose the following 95% c.l. 
ranges: 

a) -0.53 X 10"^° < (5a^ and b) 5a^ < 44.7 x 10"^°, e+e--based; (31) 
a) - 13.5 X 10"^° < 6a^ and b) 5a^ < 28.4 x 10"^°, r-based. (32) 



A lower bound on rriLSP can be derived from Eq. (31b) [Eq. (32a)] for ^ > [^ < 0]. 
Although the ^ < case can not be excluded [55], it is considered as quite disfavored 
[56] because of the poor r-decay data. 
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4.4 Collider Bounds 

For our analysis, the only relevant collider bound is the 95% c.l. LEP bound [22] 
on the lightest CP-even neutral Higgs boson, h, mass 

mh > 114.4 GeV, (33) 

which gives a [almost always the absolute] lower bound on rriLSP for ^ < [^ > 0]. 
The SUSY contributions to rrih are calculated at two loops by employing the program 
FeynHiggsFast [57] contained in the micrOMEGAs package [40]. 

5 Restrictions on the SUSY Parameters 

Applying the cosmological and phenomenological requirements given in Sec. 4, we 
delineate on the w-lsp — Afj plane the allowed parameter space of our models in 
Sees. 5.1 and 5.2 for /x > and /j, < respectively. For simplicity, as{Mz) is 
fixed to its central experimental value (equal to 0.1185) throughout our calculation. 
Note that allowing it to vary in its 95% c.l. experimental range, 0.1145 — 0.1225, 
the range of rrn,{Mz) in Eq. (1) would be slightly widened. This would lead to a 
sizeable enlargement of the allowed area for /x > due to the sensitivity of f^LSP^^ 
to the 6-quark mass (see Sec. 5.1). On the contrary, our results for /i < would 
be essentially unaffected, since the LSP annihilation to hb via an A-boson exchange 
gives a subdominant contribution to f^LSP^^ (see Sec. 5.2). 

5.1 The /i > Case 

The restrictions on the rriLSP ~ •^f2 plane with rai,{Mz) = 2.888 GeV are shown in 
Fig. 5 as solid lines, while the upper bound on rriLSP from Eq. (29) for mh{Mz) = 
2.684 GeV [mb{Mz) = 3.092 GeV] is depicted by a dashed [dotted] hue. Needless to 
say that the constraints from Eqs. (30b) and (32a) do not restrict the parameters, 
since they are always satisfied for ^u > 0. 
We observe the following: 

• The lower bounds on m-LSP are not so sensitive to the variations of mb{Mz)- 

• The lower bound on rriLSP from Eq. (33) overshadows all others. 

• The LSP annihilation via the s-channel exchange of an A-boson is by far the 
dominant (co)annihilation process near the almost vertical part of the line 
corresponding to the upper bound on w-lsp from Eq. (29). This can be 
explained by considering Fig. 6, where we draw tua and MsusY versus mLSP 
for various Afj's and the central value of mh{Mz)- We see that ula is always 
smaller than 2mLSP but close to it. We also observe that, as rriLSP or A^j 
increases, we move away from the ^-pole which, thus, becomes less efficient. 
As a consequence, J^lsp^^ increases with w-lsp or /S.f,^. 
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Figure 5: Restrictions on the thlsp — Afj plane for /i > 0. From left to right, the solid lines depict 
the lower bounds on ttilsp from Eq. (31h), (30a), (33) and the upper bound on rriLSP from Eq. (29) 
for mb{Mz) = 2.888 GeV. The dashed [dotted] line depicts the bound on rriLSP from Eq. (29) for 
mb(Mz) = 2.684 GeV [3.092 GeV]. The allowed area for mt{Mz) = 2.888 GeV is shaded. 

• The upper bound on ttt-lsp from Eq. (29) is extremely sensitive to the vari- 
ations of mb{Mz)- Especiahy sensitive is the almost vertical part of the line 
corresponding to this bound, where, as we saw above, the LSP annihilation via 
an ^-boson exchange in the s-channel is by far the dominant process. This ex- 
treme sensitivity can be understood from Fig. 7, where niA is depicted versus 
w-LSP for the lower, upper and central values of mh{Mz) in Eq. (1). We see 
that, as mi,{Mz) decreases, rriA increases and approaches 2?TiLSP- The ^-pole 
annihilation is then enhanced and f^LSP^^ is drastically reduced causing an 
increase of the upper bound on tolsp- 

• For Af2 < 0.25, bino-stau coannihilations [46] take over leading to a very 
pronounced reduction of f^LSP^^) thereby increasing the upper limit on w-lsp- 

For mb{Mz) = 2.888 GeV, the allowed ranges of rriLSP and Afj are 

176 GeV < mLSP < 615 GeV and < A^, < 1.8. (34) 



5.2 The fi <0 Case 

As is seen from Fig. 4, 2?n,LSP ^ itt-a in the /x < case. So, the LSP annihilation 
to bb via the s-channel exchange of an ^-boson is not enhanced as for /i > (the 
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Figure 6: The mass parameters tua and Msusy versus mLSP for fi > 0, mb{Mz) — 2.888 GcV 
and various values of A^j , which are indicated on the curves. 
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Figure 7: The mass parameters niA and Msusy versus m-LSP for /j. > 0, Af^ = 1 i^nd with 
mh{Mz) = 2.684 GeV (dashed lines), 3.092 GeV (dotted lines) or 2.888 GcV (solid lines). 
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Figure 8: Restrictions on the tolsp — Arj plane for /i < and mb{Mz) in the range of Eq. (1). 
The solid [dot-dashed] line corresponds to the lower bound on ttilsp from Eq. (32a) [Eq. (30h)]. 
The dashed [dotted] line corresponds to the upper bound on m-LSP [^^2] from Eq. (29). 



important annihilation channels, for /i < 0, are not only the ones with fermions // 
in the final state, but also with HZ, W^H^ and hA [43, 58]). As a consequence, the 
only available mechanism for reducing the ilLSp/i^ is the x ~ ^2 coannihilation [46] 
(for an updated study, see also Ref. [59]), which becomes efficient when Afj < 0.25. 
The restrictions from all the requirements imposed on the w,lsp — ^f2 plane 
for any rrn,{Mz) in Eq. (1) are presented in Fig. 8. The lower bound on m-LSP 
from Eq. (32a) corresponds to mb{Mz) — 3.092 GeV and is represented by a solid 
line. The maximal Af2(~ 0.032) on this line yields JIlsf^^ — 0.13 for the same 
value of mb{Mz)- As mh{Mz) decreases, the maximal Afj from Eq. (29) increases 
along the dotted line and reaches its overall maximal value at Afj — 0.035 for 
rrH){Mz) — 2.684 GeV. The upper bound on w-lsp from Eq. (29) is achieved for 
rrn,{Mz) — 2.684 GeV and corresponds to the dashed line. In the lightly shaded 
region allowed by Eqs. (32a) and (29), we find 



348 GeV < rriLSP < 680 GeV and < A^^ < 0.035. 



(35) 



The maximal m-LSP is achieved at Afj = yielding BR(6 -^ S7) ~ 5.8 x 10~^. 

On the other hand, the lower bound on rriLSP (dot-dashed line) from Eq. (30b) 
corresponds to mi,{Mz) — 3.092 GeV. In the corresponding allowed (dark shaded) 
area in Fig. 8, 

m-LSP > 1305.04 GeV (36) 
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with the minimal ttt-lsp achieved at A^j = and yielding 17lsp^ — 0.65. Note 
that the NLO corrections to BR(6 -^ 57) (see Sec. 4.2) drastically reduce this 
lower bound on m-LSP- It is worth mentioning that, in the second paper of Ref. [32] 
(which adopts the opposite sign convention for n), where the tan/3 enhanced and 
NLO SUSY QCD corrections were not included, the reduction was considerably 
higher and yielded a much less stringent restriction. 

Combining Eqs. (35) and (36), it is obvious that we are left with no simulta- 
neously allowed region. Needless to say that the constraints from Eqs. (30a) and 
(32b) do not restrict the parameters, since they are always satisfied for ^ < 0. The 
same is also valid for the bound on the lightest Higgs boson mass, Eq. (33), due to 
the heavy SUSY spectrum. 

6 The Deviation from YU 

The deviation from YU is estimated by employing Eq. (27). In the allowed (shaded) 
area of Fig. 5 (for ^ > 0) which corresponds to the central value of mh{Mz) in 
Eq. (1), the ranges of the parameters c, Sh-,-, Sh^ and tan/3 are 

0.14 < c < 0.17, 0.24 < 6hr = -5hb < 0.29 and 58 < tan/3 < 59. 

However, letting nH){Mz) vary in its 95% c.l. range in Eq. (1), we find that, in the 
corresponding allowed area, these parameters range as follows: 

0-11 < c < 0.19, 0.2 < 5hr = -6hb < 0.32 and 57 < tan/3 < 60. 

In the lightly shaded area of Fig. 8 (for /U < 0), which is obtained by allowing 
mb{Mz) to vary in the range of Eq. (1), the parameters above range as follows: 

0.01 < -c < 0.06, 0.04 < -6hr < 0.23, Sh = and 46 < tan/3 < 49. 

For rab{Mz) — 2.888 GeV, in the corresponding allowed area, these ranges become 

0.036 < -c < 0.045, 0.14 < -5hr < 0.17, 6hb = and 47 < tan/3 < 48. 

We observe that, as we increase [decrease] mb{Mz) for ^u > [/x < 0], the 
parameter |c| decreases and we get closer to exact YU. This behavior is certainly 
consistent with the fact (see also Sec. 1) that the value of mb{Mz) which corresponds 
to exact YU lies well above [a little below] its 95% c.l. range for ^ > [fi < 0]. 

Note, finally, that, for /i > 0, the required deviation from YU is not so small. In 
spite of this, the restrictions from YU are not completely lost but only somewhat 
weakened. In particular, tan/3 remains large and close to 60. Actually, our model is 
much closer to YU than generic models where the Yukawa coupling constants can 
differ even by orders of magnitude. Also, the deviation from YU is generated here 
in a natural, systematic, controlled and well-motivated way. 
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7 Direct Detection of Neutralinos 

As we showed in Sec. 5.1, our /i > model possesses a wide and natural range of 
parameters allowed by all the relevant astrophysical and experimental constraints. 
It would be, thus, interesting to investigate whether the predicted LSPs in the 
universe could be detected in the current or planned experiments. This can be 
done by first calculating the elastic scattering of the LSPs with nuclei [44, 60]. To 
accomplish this goal, we need an effective Lagrangian (see Sec. A. 2) derived from 
the MSSM Lagrangian (see Sec. A.l) and providing a reliable description of the 
neutralino— quark interaction. We also need a procedure for going from the quark 
to the nucleon level (see Sec. 7.1) and from the nucleon to the nucleus (see Sees. 7.2, 
A.3 and A.4). 

7.1 Scalar Neutralino— Proton Cross Section 

The quantity which is being conventionally used in the recent literature (see e.g. 
Refs. [61, 62, 63]) for comparing experimental [64, 65, 66] and theoretical results 
is the spin independent (SI) neutralino— proton (% — p) cross section (at zero mo- 
mentum transfer) (t|p calculated by applying Eq. (A. 24) with tIn = -^N = 1- The 
SI effective x ~ P coupling fp is calculated using the full one-loop treatment of 
Refs. [67, 68] (already used in Ref. [16]), which, however, agrees with the tree- level 
approximation (see Eq. (A. 25)) for the values of the SUSY parameters encountered 
in our model. The values in Eq. (A. 26) [62] were adopted for the renormalization- 
invariant functions P (with q = u,d, s) needed for the calculation of fp. 

Combining the sensitivities of the recent [64] and planned [65] experiments, we 
obtain the following observationally interesting region: 

a) 3 X 10"^ pb < afp and b) afp<2x 10"^ pb, (37) 

for 150 GeV < rrii^sp ^ 500 GeV. For the allowed values of the SUSY parameters of 
our model, a^ lies beyond the preferred range of DAMA [66], (1 — 10) x 10~^ pb, 
which though has mostly been excluded by other collaborations (e.g. EDELWEISS 
and ZEPLIN I [64]). 

Allowing r^LSP^^ and the hadronic inputs P to vary within their ranges in 

Eqs. (28) and (A. 26) respectively, we derive the shaded band on the ttilsp — o"|p 
plane (Fig. 9). The bold solid line corresponds to the central values of Qi.sph'^ and 
fP . Note that the width of the band is almost exclusively due to the variation of 

/^ since, for fixed ttt-lsP) (^tp is almost insensitive to the variation of f^LSP^^ within 
the range of Eq. (28) (or, equivalently, to the required variation of Afj). We observe 
that, for the allowed to-lsp's, there are o"|p's which lie within the margin of Eq. (37). 
Therefore, the LSPs predicted by our model can be detectable in the near future 
experiments. In particular, employing central values for f^LSP^^ and P , we find 

5 X 10"^ pb > 4p ^ 3 X IQ-^ pb for 179.7 GeV < rriLSP < 300 GeV. (38) 
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Figure 9: The SI x - P cross section a^p versus niLSP for p > and mh{Mz) = 2.888 GeV. 
The bold sohd hne is derived by fixing Qi^sph^ and f!^ to their central values in Eqs. (28) and 

(A. 26) respectively, whereas the shaded band by allowing Qi^sph^ and f^ to vary in their ranges 
in these equations. The region ofEq. (37), which is preferred by the various experimental projects, 
is approximately limited between the two thin solid lines. 

These values are somewhat higher than the ones obtained in similar estimations (see 
e.g. the third paper in Ref. [34]) with lower tan/3's. Furthermore, the upper bound 
on rriLSP in Eq. (34) implies a lower bound on a^ . Namely, 



afp > 1.6 (0.8) X 10 



-10 



pb. 



(39) 



where the bound in parenthesis is derived by allowing P to vary. 



7.2 Detection Rate of Neutralinos 

The total detection rate (events per day) of LSPs per kgr of detector material 
(consisting of nuclei N) can be found by [68, 69] 



^x(N) 



4 



SI p2/^ N I ^SD 'S'(Qr) 



T{Qr) 



(40) 



Here Q^ is the (recoil) energy transferred to the nucleus, Qt is the detector threshold 
energy below which the detector is insensitive to LSP— nucleus recoils, /Uj^n is the 
reduced mass (evaluated in Eq. (A. 24), a-^ is the SI [spin dependent (SD)] 
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Figure 10: The LSP detection rate RxC^Ge) for a ^^Ge detector versus jtilsp for /i > and 
mt{Mz) = 2.888 GeV. The solid line is derived by Bxing ^Lsph'^, /^'"' and Af '"' to their central 
values in Eqs. (28), (A.26) [(A.27)] and (A. 30) [(A.29)] respectively, whereas the shaded band by 



ip["] 



allowing Qi^sph , f^ and A, to vary in their ranges in these equations. 



X — N cross section (evaluated in Eq. (A. 24) [Eq. (A. 28)1) at zero Qr, p% is the total 
local CDM density and vq is the circular velocity of the sun w.r.t. the galactic rest 
frame. Given the significant uncertainties involved in the determination of the two 
last quantities, we adopt their most popular values 0.3 GeV/cm'^ and 220 km/sec 
respectively. The derivation of the SI and SD form factors F{Q^) and S'(Qr)/<S'(0) is 
also subject to considerable uncertainty. They are determined by closely following 
the analysis of Ref. [69] (and Ref. [68]). There, we can also find the function T{Q^) 
which depends on the minimal velocity of the incident LSPs and the velocity of 
earth (w.r.t. the galactic rest frame), which generates an annual modulation [70]. 

In Fig. 10, we display R^i^^Ge) as a function ttt-lsp for Qt = H keV and for the 
second of June. The shaded region is derived by allowing r^Lsp/i^, pN and ISlg 
(see Sees. A. 3 and A. 4) to vary within their ranges in Eqs. (28), (A.26) [(A.27)] and 
(A. 30) [(A.29)] respectively, whereas the solid line by fixing them to their central 
values. Note that the width of the band is almost exclusively due to the variation of 
(t|J,j, since (t|^ is almost insensitive to the aforementioned variations and an order 
of magnitude smaller than (t|^. The resulting Ry^(J^Ge) is a little lower than the 
one obtained with the naive choice F{Qr) = S{Qr)/S{0) = 1 and significantly lower 
than the one corresponding to a wino or Higgsino LSP (compare with Ref. [71]). 
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8 A Resolution of the fi Problem 

An important shortcoming of MSSM is that there is no understanding of how the 
SUSY |U term, with the right magnitude of |/i| ~ 10^ — 10^ GeV, arises. One way 
[23] to solve this fi problem is via a PQ symmetry U{1)pq [24], which also solves the 
strong CP problem. This solution is based on the observation [72] that the axion 
decay constant fa, which is the symmetry breaking scale of U{1)pq, is (normally) 
intermediate (~ 10^"*^ — 10"*^^ GeV) and, thus, |/i| ~ f^/Ms. The scale fa is, in turn, 
~ (^^3/2^3)^) where 7713/2 ~ 1 TeV is the gravity-mediated soft SUSY breaking 
scale (gravitino mass). In order to implement this solution of the /i problem, we 
introduce [23] a pair of superfields N and N (see Table 1) with the following non- 
renormalizable couplings in the superpotential [73]: 

Here, A^, X'^, A and Apq are taken positive by redefining the phases of N and N. 
After SUSY breaking, the N'^N'^ term leads to the scalar potential 

I [(|iV|-|iV|)2 + 2|Ar||iV|] 



^PQ = ^2/2 + 4Ai>Q 



NN 



Ms 



liViVp 
+2|A|m3/2ApQ^^cos(e + 2^ + 20), (42) 

where A is the dimensionless coefficient of the soft SUSY breaking term corre- 
sponding to the superpotential term N'^N'^ and e, 9, 6 are the phases of A, N , N 
respectively. Minimization of Vpq then requires |A^| = \N\, e + 26 + 29 = tt and Vpq 
takes the form 

For 1^1 > 4, the absolute minimum of the potential is at 



fa _ I ^ |A| + VI^P-12 



KiV)| = |(Ar)| ^ :^ = ^m3/2Ms y ^ ^m^/2Ms. (44) 

The // term is generated predominantly via the terms N'^h? and N'^h''^ of Eq. (41) 
with |/i| ~ \{N)\'^ /Ms, which is of the right magnitude. 

The potential Vpq also has a local minimum at A^ = iV = 0, which is separated 
from the global PQ minimum by a sizable potential barrier preventing a successful 
transition from the trivial to the PQ vacuum. This situation persists at all cosmic 
temperatures after reheating, as has been shown [27] by considering the one-loop 
temperature corrections [74] to the scalar potential. We are, thus, obliged to assume 
that, after the termination of inflation, the system emerges with the appropriate 
combination of initial conditions so that it is led [75] to the PQ vacuum. 
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9 Proton Stability 

We can assign baryon number 1/3 [—1/3] to all color triplets [antitriplets] . Recall 
that there are (anti)triplets not only in F, F'^ but also in H^, H^, G. Lepton number 
is then defined via B — L. Before the inclusion of the extra Higgs superfields (see 
Table 1), B (and L) violation comes from the following terms [27] 

pcpcjjcjjc^ FFH^H^hh, FFH'^H'^N'^ (45) 

(and terms containing the combinations (H'^)'^, (H'^)'^), which give couplings like 
u'^(F(F^v'fj (or u^d^UffCff), uddfjjv'fj (or udu'jje'j^) with appropriate coefficients. Also, 
the terms GH'^H'^ and GH'^H^ give rise to the B (and L) violating couplings 
g'^u'jjd'^, g'^u'jjd'^jj. All other combinations are B (and L) conserving since all their 
SU{4)c 4's are contracted with 4's. 

The dominant contribution to proton decay comes from effective dimension five 
operators generated by one- loop diagrams with two of the n^, d'jj or one of the n^, 
d'jj and one of the i^|^, e^ circulating in the loop. The amplitudes corresponding to 
these operators are estimated to be at most of order 7713/2 Mgut/A^s ^ 10~^^ GeV^^. 
This makes the proton practically stable. 

After the inclusion of the superfields h' and h' , the couplings 

FFH^H^hh', FFH^H''h'h' (46) 

(as well as the new couplings containing arbitrary powers of the combinations {H'^) , 
{H^)^) give rise [9] to additional B and L number violation. However, their con- 
tribution to proton decay is subdominant to the one arising from the terms of Eq. 
(45). One can further show [9] that the inclusion of the superfields (j) and (f) also 
gives a subdominant contribution to the proton decay. 

10 The Inflationary Scenario 

One of the most promising inflationary scenarios is hybrid inflation [76] , which uses 
two real scalars: one which provides the vacuum energy density for inflation and a 
second which is the slowly varying field during inflation. This scheme is naturally 
incorporated [77] in SUSY GUTs (for an updated review, see Ref. [78]), but in its 
standard realization has the following problem [79]: if the GUT gauge symmetry 
breaking predicts monopoles (and this is the case of Gps which predicts doubly 
charged monopoles [80]), they are copiously produced at the end of inflation leading 
to a cosmological catastrophe [81]. One way to remedy this is to generate a shifted 
inflationary trajectory so that Gps is already broken during inflation. This could 
be achieved [27] in our SUSY GUT model even before the introduction of the extra 
Higgs superfields (see Table 1), but only by utilizing non-renormalizable terms. The 
inclusion of h' and h' does not change this situation, which thus also holds in our 
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/i < model. This model is though excluded by imposing the restrictions of Sec. 4. 
On the other hand, for n > 0, the inclusion of 4> and (f> very naturally gives rise 
[25] to a shifted path, but only with renormalizable interactions (for similar recent 
analyses in the context of a SU{5) SUSY GUT, see Ref. [82]). 

10.1 The Shifted Inflationary Path 

As we showed in Sec. 2.2, the deviation from YU which is needed for ^ > can be 
achieved by the inclusion of the superfields (p, (/). In the presence of these superfields, 
a new version [25] of shifted hybrid inflation [27] can take place, without invoking 
any non-renormalizable superpotential terms. 

Indeed, these fields lead to three new renormalizable terms in the part of the 
superpotential which is relevant for inflation. This is given by 

W = kS{H^H'' - M^) - pScj)"^ + m<j)4) + X^H^H^, (47) 

where M,m ~ MquT; and k, (3 and A are dimensionless coupling constants with 
M, m, K, A > by field redefinitions. For simplicity, we take /? > (the parameters 
are normalized so that they correspond to the couplings between the SM singlet 
components of the superfields). 

The scalar potential obtained from W is given by 

V = \k{H^H^ - M^) -(3(l)^f + \2(3Sct) - m^^ + |m0 + XH^H^f 

+ \kS + X4)\^ {\H^\^ + |F"|2) + D - terms. (48) 

Vanishing of the D-terms yields H'^* = e^'^H'^ {H'^, H'^ lie in the v'jj, v'fj direction). 
We restrict ourselves to the direction with i? = which contains the shifted infla- 
tionary path and the SUSY vacua (see below). Performing appropriate R and gauge 
transformations, we bring H'^, H^ and S to the positive real axis. 

From the potential in Eq. (48), we find that the SUSY vacuum lies at 

^^(5r-^(-v^)'i-/l(5r («> 

with S* = and (j) = 0, where ^ = [5}? M'^ / Km? < 1/4. The potential possesses a 
'shifted' flat direction (besides the trivial one) at 



M2 \mJ 2(At2 + A2) ' M 2V/3e A 

with S > and a constant potential energy density Vq given by 



M4 k2 + A2 V4^ 
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which can be used as inflationary path. Vq 7^ breaks SUSY on this path, while 
the constant non-zero values of H'^, H'^ break the GUT gauge symmetry too. The 
SUSY breaking implies the existence of one-loop radiative corrections [83] which lift 
the classical flatness of this path yielding the necessary inclination for driving the 
inflaton towards the SUSY vacuum. 

The one- loop radiative corrections to V along the shifted inflationary trajectory 
are calculated by using the Coleman- Weinberg formula [84]: 

i 

where the sum extends over all helicity states i, F-i and Mf are the fermion number 
and mass squared of the ith state and A is a renormalization mass scale. In order 
to use this formula for creating a logarithmic slope which drives the canonically 
normalized real inflaton field a = \J2(^K^ + A^)5'/A towards the minimum, one has 
first to derive the mass spectrum of the model on the shifted inflationary path. This 
is a quite complicated task and we will skip it here. 

10.2 Inflationary Observables 

The slow roll parameters are given by (see e.g. Ref. [85]) 

..^(my„a,.,4^, (.3, 

where the primes denote derivation w.r.t. the real normalized inflaton field a and 
7TT-P ~ 2.44 X 10^^ GeV is the reduced Planck scale. The conditions for inflation to 
take place are e < 1 and |?7| < 1. 

The number of e-foldings that our present horizon scale suffered during inflation 
can be calculated as follows (see e.g. Ref. [85]): 






^Q^^J, ^^^ ^ 1" ( 4-41 X 10^^ T/ V,^ ] , (54) 



where af [ctq] is the value of a at the end of inflation [when our present horizon scale 
crossed outside the inflationary horizon] and T^ ~ 10^ GeV is the reheat temperature 
taken to saturate the gravitino constraint [86]. 

The quadrupole anisotropy of the CMBR can be calculated as follows (see e.g. 
Ref. [85]): 

'^T\ ^ 1 VW 



TJq 12V5mlV'{aQ) 



(55) 



Fixing {6T/T)q ~ 6.6 x 10 ^ to its central value from COBE [26] (under the con- 
dition that the spectral index n = 1), we can determine one of the free parameters 
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(say (3) in terms of the others (m, k and A). For instance, one finds [25] /3 = 0.1 
for m = 4.35 x 10^^ GeV and k = A = 3 x 10^^. In this case, the instability- 
point of the shifted path is at ac ~ 3.55 x 10^'^ GeV, ct/ ~ 1.7 x lO^'^ GeV and 
aq ~ 1.6 X lO^s GeV. Also, M ~ 2.66 x 10^^ GeV, Nq ~ 57.7 and the spectral 
index n ~ 0.98. Note that the slow roll conditions are violated and inflation ends 
well before reaching the instability point at ac- We see that the COBE constraint 
can be easily satisfied with natural values of the parameters. Moreover, superheavy 
SM non-singlets with masses <C MquTi which could disturb the unification of the 
MSSM gauge coupling constants, are not encountered. 

10.3 SUGRA Corrections 

As we emphasized, the new shifted hybrid inflation occurs at values of a which 
are quite close to the reduced Planck scale. Thus, one cannot ignore the SUGRA 
corrections to the scalar potential. 

The scalar potential in SUGRA, without the D-terms, is given by 



Y ^ ^K/ml 






(56) 



where K is the Kahler potential, Fi = Wi + KiW/nip, a subscript i [i*] denotes 
derivation w.r.t. the complex scalar field </>* [</>* *] and K'' ^ is the inverse of the 
matrix Kji*. 

Consider a (complex) inflaton S corresponding to a flat direction of global SUSY 
with Wjs = 0. We assume that the potential on this path depends only on |S|, 
which holds in our model due to a global symmetry. From Eq. (56), we find that 
the SUGRA corrections lift the flatness of the S direction by generating a mass 
squared for S (see e.g. Ref. [87]) 

-'- = ;^ - ^ + Ei^^r^g^w, + . . , (57) 



mf, TJip 



«j 



where the right hand side (RHS) is evaluated on the flat direction with the explicitly 
displayed terms taken at S = 0. The ellipsis represents higher order terms which 
are suppressed by powers of |S|/?7ip. The slow roll parameter ij then becomes 

^ = 1 - ^ + 9^ Y.m*Kg^, + . . . , (58) 

which, in general, could be of order unity and, thus, invalidate [77, 88] inflation. 
This is the well-known rj problem of inflation in local SUSY. Several proposals have 
been made to overcome this difficulty (for a review, see e.g. Ref. [87]). 

In all versions of SUSY hybrid inflation, there is an automatic mutual cancella- 
tion between the first two terms in the RHS of Eq. (58). This is due to the fact that 
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Wn = on the inflationary path for all field directions n which are perpendicular to 
this path, which implies that |VFeP = Vq on the path. This is an important feature 
since, in general, the sum of the first two terms in the RHS of Eq. (58) is positive 
and of order unity, thereby ruining inflation. Moreover, Wn = also implies that 
the only contribution to the third term in the RHS of Eq. (58) comes from the term 
in K which is quartic in S. So the third term can be suppressed by mildly tuning 
just one parameter [89] and inflation could remain intact provided that the terms 
in the ellipsis can be ignored. 

However, in our present model, inflation takes place at values of |S| close to nip. 
So, the terms in the ellipsis in the RHS of Eq. (58) cannot be ignored (in contrast to 
the case of the old version [27] of shifted hybrid inflation) and may easily invalidate 
inflation. We, thus, need to invoke here a mechanism which can ensure that the 
SUGRA corrections do not lift the flatness of the inflationary path to all orders. A 
suitable scheme has been suggested in Ref. [90]. It has been argued that special 
forms of the Kahler potential can lead to the cancellation of the SUGRA corrections 
which spoil slow roll inflation to all orders. In particular, a specific form of i^(S) 
(used in no-scale SUGRA models) was employed and a gauge singlet field Z with a 
similar K{Z) was introduced. It was pointed out that, by assuming a superheavy 
VEV for the Z field through D-terms, an exact cancellation of the inflaton mass on 
the inflationary trajectory can be achieved. 

The form of the Kahler potential for S used in Ref. [90] is given by 

i^(|S|2) = -AAmi>ln(l--i^), (59) 

where AA = 1 or 2; here we will take N = 2. In this case, the kinetic term of the 
real normalized inflaton field a (note that |S| = cr/^/2) is {l/2){d'^K/dT.dT.*)a'^ , 
where the overdot denotes derivation w.r.t. the cosmic time t and d'^K/d'EdT,* = 
(1 — cr^/2AA?72p)~^. Thus, the Lagrangian density on the shifted path is given by 



C = a^{t) 



^■21- cr 



2 \-2 



o'^M 1 - TTTTzrr -n^ 



2M' 



m 



(60) 



where a{t) is the scale factor of the universe. 

The evolution equation of a is found by varying this Lagrangian w.r.t. a: 



d + 3iJd + aM 1 



•2 



^2 X a 



2J\fmf. J J\fmf> 



a^ 



l-TTT^) +V\ct)=0, (61) 
P- 



2Mm^ 



where H is the Hubble parameter. During inflation, the 'friction' term 3Ha dom- 
inates over the other two terms in the brackets in Eq. (61). Thus, this equation 
reduces to the modified inflationary equation 
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Note that, for a <^ \/2Mmp , this equation reduces to the standard inflationary 
equation. 

To derive the slow roll conditions, we evaluate the sum of the first and the third 
term in the brackets in Eq. (61) by using Eq. (62): 



, + ,M:_^V"^ = IV)«V)^A ^' 



2MmlJ Mml 3H^ \ 2Mml 

^ 1 - TTTTzrr + ^77^^ 1 - TTTTzrr TTzrr • (63) 



?,H " V^ 2J\fmlJ 3H \ 2Mml) Mml 

Comparing the first two terms in the RHS of Eq. (63) with Ha, we obtain 

1 2 fV'ia)V f. a^ ^^ 



'-r^A-^) I'-w;^' -'' ^''^ 



\r]\ ~ m 



V"{cy\ 



Vr 







a^ 



The third term in the RHS of Eq. (63), compared to Ha, yields \/2ae^''^ /Nm-p < 1, 
which is automatically satisfied provided that Eq. (64) holds and a < Mra-pl\f2. 
The latter is true for the values of a which are relevant here. We see that the slow 
roll parameters e and r\ now carry an extra factor (1 — a"^ /2J\fm'^)'^ < 1. This leads, 
in general, to smaller cr/'s. However, in our case, u/ <C v2J\fmp (for J\f = 2) and, 
thus, this factor is practically equal to unity. Consequently, its infiuence on af is 
negligible. 

The formulas for Nq and {5T/T)q are now also modified due to the pres- 
ence of the extra factor (1 — a'^ /2J\frn^)'^ in Eq. (62). In particular, a factor 
(1 — o"^/2AAmp)~^ must be included in the integrand in the RHS of Eq. (54) and a 
factor (1 — crQ/2AAmp)~'^ in the RHS of Eq. (55). We find that, for the a's under con- 
sideration, these modifications have only a small infiuence on aq if we use the same 
input values for the free parameters as in the global SUSY case. On the contrary, 
{5T /T)q increases considerably. However, we can easily readjust the parameters so 
that the COBE requirements are again met. For instance, {5T/T)q ~ 6.6 x 10^® 
is now obtained [25] with m = 3.8 x 10^^ GeV keeping k = A = 3 x 10"^, /3 = 0.1 
as in global SUSY. In this case, ac ^ 2.7 x 10^^ GeV, af ~ 1.8 x lO^^^ GeV and 
aq ~ 1.6 X 10^^ GeV. Also, M ~ 2.6 x 10^^ GeV, Nq ~ 57.5 and n ~ 0.99. 

11 Conclusions 

We studied the CMSSM with Aq = applying a suitable set of YQUCs which 
originate from SUSY GUT models based on Gps. For each sign of ^, an appropriate 
YQUC was chosen from this set so that an adequate deviation from YU which allows 
an acceptable mb{Mz) is ensured. We, also, imposed the constraints from the CDM 
in the universe, b — > 57, 5a^ and rrih- We concluded that 
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• For ;U > 0, there exists a wide and natural range of CMSSM parameters 
which is consistent with all the above constraints. We found that tan f3 ranges 
between about 58 and 61 and the asymptotic splitting between the bottom (or 
tau) and the top Yukawa coupling constants varies in the range 26 — 35% for 
central values of mi,{Mz) and as{Mz)- 

• For /U < 0, despite the fact that, considering the r-based calculation for af^, 
the (5q^ and CDM criteria can be simultaneously satisfied, the model is ex- 
cluded since the b ^ sj and CDM requirements remain incompatible. How- 
ever, the deviation from YU needed for correcting mi,{Mz) is much smaller 
than in the previous case. 

The predicted LSP mass in the fi > model can be as low as about 176 GeV 
with a^ in the range of the sensitivity of the planned direct CDM detectors, though 
with a low detection rate. Moreover, the fi > model resolves the /i problem of 
MSSM, predicts stable proton and gives rise to a new version of shifted hybrid 
inflation. The inflationary scenario relies on renonormalizable terms only, can be 
consistent with the COBE constraint on the CMBR with natural values for the 
relevant parameters and avoids overproduction of monopoles at the end of inflation. 
Inflation takes place along a classically flat direction, where Gps is spontaneously 
broken to Gsm- A readjustment of the input values of the free parameters is required 
after the introduction of a specific Kahler potential and an extra gauge singlet with 
a superheavy VEV via D-terms, which are to be included in order for the SUGRA 
corrections not to invalidate inflation. 
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Appendix A 

Neutralino— Nucleus Elastic Cross Section 

In this Appendix, we sketch the derivation of the tree-level approximation iT?p^ to 
the total cross section at zero momentum transfer for the elastic scattering process 
xN ^ xN with N being a nucleus target. This quantity can be split into a scalar 
7?}^) and an axial- vector or SD ((t|^) 



or SI (o"|}^) and an axial- vector or SD ((t|^) part [60, 61, 62, 67, 



^0 — ^si , ^SD /' A 1 ^ 



In the following, we present the ingredients needed for the calculation of each part. 
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A.l The Relevant MSSM Lagrangian 

We present, in terms of mass eigenstates, all the Feynman rules which are necessary 
for studying the neutralino-quark elastic scattering process. 

A. 1.1 Neutralino—Squark— Quark Vertex (Ref. [91], Fig. 24) 

The Lagrangian which describes the neutralino—squark— quark interaction is 

2 

C^gg = V2gtwqy( qIlPl + qIr Pr ) m + h.c. , (A.2) 



where h.c. denotes the hermitian conjugate, Pl[r] = (1 ~ [+]75)/2, <? is a quark field, 
g is the SU(2)l gauge coupling constant and tw = t&nOy/. Here, we considered 
general flavor-diagonal squark mixing and gi, q2 are the squark mass eigenstates 
defined as follows: 

(I) --M£)'*-^-(:: ?]• (^^^> 

Sq = sin 9q and Cq = cos 6q with 9q being the ql — Qr mixing angle. The coefficients 
g^^ and gfj^ in the mass eigenstate basis are 

dti = Cg gli + «9 91r, 9\r = Cq g^L + Sg 9rr^ (A-4) 

92L = ~'^<? 9lL + ^Q 9lRi 92R = ~-^q 9rl + ^q Qrr ■ (A. 5) 

The gill Blr^ 9rr ^^d g'j^j^ coefficients for the up-type quarks n, c and t are 

and for the down-type quarks d, s and b are 

where svv' = sin 6*14^, cw = cos ^vi^, s^ = sin/3, c^ = cos/9, m„[^ are the masses of the 
up-type [down-type] quarks and N is the matrix which diagonalizes the neutralino 
mass matrix. Finally, N[-^ and N[2 are defined as follows [91]: 
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A.1.2 Z-boson-Quark-Quark Vertex (Ref. [92], Fig. 71) 

The Lagrangian which describes the Z-boson— quark— quark interaction is 

Czqg = gzq7^[RgPR + LgPL)qZ^, (A.ll) 

where Rq and Lq for the up-type [down-type] quarks read 

Luldi = -[+] U - [+] 2 ^-y^^H^) and R^idi = -2 '^^^^w, (A.12) 

and gz = g/2cw- 

A.1.3 Z-boson— Neutralino-Neutralino Vertex (Ref. [92], Fig. 75) 

The Lagrangian which describes the Z-boson— neutrahno—neutrahno interaction is 

^Zxx = 9ZxxX7''l5XZf,, where fif^^^ = Y ('^i^'^~ '^^'''^) • (^-13) 

A. 1.4 Higgs-Quark-Quark Vertex (Ref. [91], Fig. 7) 

The Lagrangian which describes the Higgs— quark— quark interaction is 

^h,H,Aqq = qyghqqh + gUqqH + igAqq^bAjq, (A.14) 

where the ghqq,gHqq and gAqq coefficients for the up- and down- type quarks are 

_ g rriuCq _ grudSg (\^r,\ 

2MwSi3 2MwCi3 

gAuu = - 2M ^' ^^'^'^ ^ ~ 2M *^^^ ' (A- 17) 

where c^ = cos a and s^ = sina with a being the Higgs mixing angle [91]. 

A. 1.5 Higgs— Neutralino-Neutralino Vertex (Ref. [91], Fig. 21) 

The Lagrangian which describes the Higgs— neutralino—neutralino interaction is 

^KH,A XX = 2 ^ \9hxx ^ + 9Hxx ^ + ^9Axx 75 ^ j X , (A. 18) 

where 

gh[H\xx= 9 [s a [-Ca\Q 11+ Ca[Sa]Sll\ and gAxx = 9\Sf5 Qll - Cp Suj (A.19) 

with 

Qii = A^i3 {Ni2 - twNii) and 5ii = TVm (^12 - twNu) . (A.20) 



32 



G. Lazarides and C. Pallis 




Qia 




(a) 





Figure 11; The five Feynman diagrams contributing to the neutralino— quark elastic scattering 
process xQ ~* XQ- (^) ^he two diagrams with squarli qi_2 exchange, (b) the two diagrams with 
spin-0 neutral-Higgs-boson {h,H) exchange, and (c) the diagram with spin-1 Z-boson exchange. 

A. 2 Neutralino— Quark Effective Lagrangian 

From the expressions above, it is possible [60, 67, 93], after the appropriate Fierz 
rearrangement, to derive the coefficients a^ and a^ contained in the effective four- 
fermion Lagrangian for the neutrahno— quark elastic scattering: 



£ 



eflF 



af ixx) (qq) + a 



SD 



{xi^i^x) {qipiim) + 



(A.21) 



where the ellipsis represents (i) terms (such as X75X9759) which generate velocity- 
dependent and, thus, negligible [93] contributions to the elastic neutralino— quark 
cross section, (ii) loop corrections of order m~^ which were previously considered 
in the literature [67] and can be safely neglected for rriq^.^ 3> m^ [61, 62] (see 
below). The five tree- level Feynman diagrams contributing to the neutralino— quark 
scattering process are shown in Fig. 11. The s-channel squark-exchange diagrams 
contribute to both a^ and a^ , while the i-channel Higgs-exchange diagrams only to 
aq. The t-channel Z-exchange diagram contributes only to a^. These coefficients, 
in the non-relativistic limit, read 



a 



SI 



ft'w E 



a 



SD 



9iL 9iR 



2 

1,2 1^ 



{my, + mq)2 



9hxx 9hqq 9HxX 9Hqq 






E 



9!l\ +1 9 



iR\ 



1,2"^! 



[rrii 



rrir, 



2ml 



2m]j 



2 9Zxx(Rq - Lq) 



2M: 



(A.22) 
(A.23) 



Finally, let us note that, besides the neutralino— quark interaction in Eq. (A.21), 
there are also neutralino— gluon interactions with contributions arising from one-loop 
heavy quark and Higgs diagrams and the so-called twist-2 operators [67, 68]. How- 
ever, these corrections are negligible for the values of SUSY parameters encountered 
in our /i > model (see Sec. 7). 
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A. 3 Neutralino— Nucleus Spin Independent Cross Section 

The SI part of a^j^ can be parameterized as 

^^In = - /"In ( -^n/p + (^N - ^n) fn ) , where ^j^n = -, (A.24) 

is the neutrahno— nucleus reduced mass and Zn and vIn denote the atomic number 
and the mass number of the nucleus respectively. In the limit of rrij^ ^ niq-^^ ^ (which 
is, in general, true unless q2 is the next-to-LSP coinciding with the lightest stop 
[47, 94] or sbottom [16] quark) and ruq <C friq^.^ , the effective coupling /p[„] of the 
LSP to proton [neutron] is given by 

/pM= E ^^/?"'«f + :|/?^ E ^^«f (A-25) 

/^L J A-,. ^ T, y 27 ^G ^-^ rn„ ^ 

q=u,d,s ^ q=c,b,t 

to lowest order in m^ , where mpin] is the proton [neutron] mass and the parameters 
P and f" are taken to be [62] 

fP^ = 0.02 ± 0.004, fP = 0.026 ± 0.005, /^^ = 0.118 ± 0.062, (A.26) 

/" = 0.014 ± 0.003, /" = 0.036 ± 0.008, /" = 0.118 ± 0.062. (A.27) 

iu J-d ^s 

with /^'"^ = 1 — X](7=M rf s /t • I^ our numerical analysis, we use the running masses 
of h and t quark at the scale rrij^. 

A. 4 Neutralino— Nucleus Spin Dependent Cross Section 



The SD part of u^i^ can be parameterized as follows: 



,SD _ 16 



^^N = — 4n A^ <^n( Jn + 1), where A = — (^ Ap {Sp) + A„ (5„)) (A.28) 



and Jn is the total angular momentum of the nucleus (9/2 for '^^Ge). Also, (Spr^i) 
is the expectation value of the spin content of the proton [neutron] group in the 
nucleus with explicit value, for a '^^Ge target, {Sp)^^ = 0.011 [{Sn).j.^ = 0.491] in 
the shell model. Finally, the coefficient Apr„i is parameterized in terms of the quark 

spin content of the proton [neutron] A^ and the effective axial-vector couplings 
<Xq^ (<? = u, d, s) as 

Kin] = E «f ^?"^ ' ^h«^« ^" = ^2' ^d = K and A^ = A? (A.29) 

q=u,d,s 

with the factors Ag taken to be [62] 

AP = +0.78 ± 0.02, A^ = -0.48 ± 0.02 and A^ = -0.15 ± 0.02. (A.30) 

Needless to say that the formalism of this section remains valid even in the case of 
light stop [47, 94] or sbottom [16] quarks. 
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